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Abstract 

We have investigated the static properties of one-dinrensional planar Josephson tunnel junctions 
in the most general case of elliptic annuli. We have analyzed the dependence of the critical current in 
the presence of an external magnetic field applied either in the junction plane or in the perpendicular 
direction. We report a detailed study of both short and long elliptic annular junctions having 
different eccentricities. For junctions having a normalized perimeter less than one the threshold 
curves are derived and computed even in the case with one trapped Josephson vortex. For longer 
junctions a numerical analysis is carried out after the derivation of the appropriate Perturbed sine- 
Gordon Equation. For a given applied field we find that a number of different phase profiles exist 
which differ according to the number of fluxon-antifluxon pairs. We demonstrate that in samples 
made by specularly symmetric electrodes a transverse magnetic field is equivalent to an in-plane 
field applied in the direction of the current flow. Varying the ellipse eccentricity we reproduce 
all known results for linear and ring-shaped Josephson tunnel junctions. Experimental data on 
high-quality Nb/Al-AlOx/Nb elliptic annular junctions support the theoretical analysis provided 
self-field effects are taken into account. 

PACS numbers: 74.50.+r, 41.20.-q, 03.75.Lm 
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I. INTRODUCTION 


Circular annular Josephson tunnel junctions (JTJs), consisting of two superconducting 
rings coupled by a thin dielectric tunneling layer were recognized to be the ideal benchmark 
to test both the statics and the dynamics of sine-Gordon solitons in the presence of a periodic 
potential 1 " 1 . In this context, a soliton is a current vortex, also called a Josephson vortex 
or a fluxon, circulating around the junction and carrying one magnetic flux quantum. A 
spatially periodic potential for the fluxon can be easily implemented by a uniform magnetic 
field applied in the plane of the annulus. However, for JTJs having a not simply-connected 
topology the most general and at the same time regular geometry is provided by an elliptic 
annulus. At variance with a circle that has infinitely many axes of symmetry, an ellipse has 
two axes of symmetry. This geometrical symmetry breaking which comes with an associated 
non-uniformity of the radius of curvature has been exploited in several physical systems, e.g., 
to increase the focal depth 5 ^ or the resolving power 6 of annular apertures and as antenna 
reflectors in the microwave region' . Elliptic annular Josephson tunnel junctions (EAJTJs) 
serve as an handy tool for the realization of complex periodic potentials, including those 
lacking spatial reflection symmetry, known as ratchet potentials-®! An additional motivation 
to study EAJTJs is the intention to cast in one unique class many apparently different JTJ 
conhgurations, including the linear geometry commonly studied in the context of JTJs. 

In this article we focus our attention on the static conhgurations of the phase in EAJTJs; 
the dynamics of solitons will be the subject of another article. Most of the work will be 
focused on EAJTJs having the so called Lyngby-type geometry^, that refers to a specularly 
symmetric configuration in which the height of the current carrying electrodes matches one 
of the ellipse outer axis (e.g. see Eigures |TT|(a) and (b)). We have measured the dependence 
of the zero-voltage Josephson current for a large number of samples. The experimental data 
are compared with the analytical results and with the numerical results obtained by solving 
the appropriate partial differential equation. 

The paper is organized as follows. In next Section we first consider a quater-clliptic annular 
junction immersed in a uniform in-plane magnetic field and derive the threshold curves 
for junctions having different ellipticity; later we extend the analysis to full ellipses with 
possible Josephson vortices trapped in the annular barrier. In Section III we derive the 
appropriate partial differential equation for an electrically long EAJTJ; later we present the 
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results of the numerical simulations and outline the effects of the magnetic held induced 
by the junction bias current. Thereafter, we investigate the consequences of a magnetic 
held applied perpendicular to the junction plane for different geometrical configurations; for 
Lyngby-type samples we establish the equivalence between a transverse magnetic held and 
an in-plane magnetic held applied along the how of the bias current. Besides, we suggest a 
simple geometrical configuration that implements the ideal step-like deterministic periodic 
ratchet potential. In Section V we describe the fabrication of our Nb/Al-AlOx/Nb samples, 
the different geometries that have been realized and the experimental setup; finally, we 
present and discuss the experimental data of long EAJTJs with both in-plane and transverse 
magnetic helds. The conclusions are drawn in Section VI. 


II. SHORT ELLIPTIC ANNULAR JUNCTIONS 

In this section we derive the equations which describe the behavior of a small EAJTJ in 
the presence of a uniform (static) magnetic held, H y, applied along one of its axis. In order 
to confer the largest generality to the analysis, we will begin by considering a quarter-elliptic 
planar junction which, depending on its eccentricity, will include also the cases of linear and 
quarter-annular^ junctions. Later on, we will treat an EAJTJ as the parallel combination 
of four quarter-elliptic ones subject to periodic boundary conditions. 


A. Quarter-elliptic junctions 

As shown in Figure [lj let the quarter-ellipse lay in the plan identified by the X-Y Carte¬ 
sian coordinate system whose origin coincides with the center of symmetry and whose axes 
are directed along the principal semi-axes a and b of its master ellipse. The curve is de¬ 
scribed by the the parametric equations x = asinr and y = b cost, where r (0 < r < 7 t/ 2) 
is a parameter measured clockwise from the positive Y'-axis such that r = ArcTan px/y , not 
to be confused with the polar angle 6 = ArcTan x/y. We defined p = b/a the axes ratio 
and we will also make use of the master-ellipse eccentricity e 2 = 1 — p 2 . For a circular arc, 
p = 1 (no eccentricity), so r and 6 coincide, while for p ^ 1, r = 9 only for 6 = mn/2. 
In the case of a thin circular ring with mean radius r, we would introduce the curvilinear 
coordinate s(9) = r6 such that s linearly increases by one perimeter (circumference) as 6 
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FIG. 1. Quarter-elliptic Josephson planar tunnel junction in a uniform in-plane magnetic field 

H|| = (0,F|j). 

changes by 2n. Along an elliptic arc we instead introduce the non-linear curvilinear coordi¬ 
nate s(r) = a fJZ(r')dr' = aE(r, e 2 ), where I(r) = \J 1 — e 2 sin 2 r = \J cos 2 r + p 2 sin 2 r is 
the integrand of the incomplete elliptic integral of the second kind, E(r, e 2 ). 

In Josephson’s original description the quantum mechanical phase difference, 0, across the 
barrier is related to the magnetic field, H, inside the barrier 11 : 

V0 = k.H x lu, (1) 

in which u- is a unit vector orthogonal to the junction plane and kG 1 = ^ 0 /2np 0 d m , where 
$o is the magnetic flux quantum, /i 0 the vacuum permeability, and d m = Xb tarih db/2Xb + 
X t t&nhd t /2\ t is the junction magnetic penetration depth®!, w here A b,t and d^t are, respec¬ 
tively, the bulk magnetic penetration depths and thicknesses of the junction bottom and 
top films; d m reduces to A& + X t in the case of thick superconducting films (d^t larger than 
~ 4A b,t). The subscripts b and t will be adopted to label the quantities which refer, respec¬ 
tively, to the bottom and top/wiring layers. 

H in Eq.Q is the total magnetic field that, in general, is given by the sum of an externally 
applied field and the self-held generated by the current flowing in the junction. Eq.Q 
states that, among other things, 0 is not sensitive to fields along the Z-direct ion; this is 
only true in the ideal case of a bare Josephson sandwich with no electrical connections. In 
real devices the presence of the current carrying electrodes not only alters the effect of a 
barrier-parallel held, but also makes the junction sensitive to a transverse held 1314 . The 
way such held induces screening currents, which, in turn, generate in-plane magnetic helds 
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is well understood in rectangular JTJs; in Section IV, we will discuss how a similar approach 
works with not simply connected junctions. 

For the time being, we assume that the quarter-elliptic junction is electrically short, 
i.e., the arc mean length, L , is small compared to the Josephson penetration length 
A j = \J < ho/27r/i 0 dj J c , where dj is the junction current thicknesd^ dj = A&cothdb/Ab + 
AiCothdt/A t > Xb + Xt > d m ; for thick film junctions, dj ~ d m ~ A& + A t . Further, we 
assume that the arc has a finite width, W, much smaller than the quarter-ellipse semi¬ 
axes. Under these conditions a spatially homogeneous field applied in the X-Y plane fully 
penetrates into the barrier. In our specific case, the external field is applied perpendicular 
to the a-axis, i.e., along the positive Y -direction, Hy = (0, ify). According to Eq. ([Tj) , the 
Josephson phase only changes in the X-direction, i.e., d(j)/dy = 0 and d(f)/dx = —kH y; then 
(p(x) = —kH ||a: + 00) where 0o is an integration constant. The dependence on r is: 

0(t) = —h\\ sin r + 0 O , (2) 


where h\\ is the strength of the in plane field, H y, normalized to (kci) - 1 = < Fo/27r/io d m a. The 
local density, Jj, of the Josephson current at a point r inside the barrier area is 11 Jj(r) = 
J c (r) sin0(r), where J c is the maximum Josephson current density. The Josephson current, 
Jj, through the barrier is obtained integrating Jj over the junction area, Ij = f A JjdA ; in 
force of the one-dimensional approximation, A = WL and the elementary surface element is 
gL4 = dWdL, where dL = ds = aX(r)dr is the elliptic elementary arc. Assuming J c uniform 
over the barrier area and recalling that 0 is constant along the annulus width: 


r r 7r / 2 

Ij = J C W /sin0(r)dL(r) = J c Wa / I(r) sin(— h\\ sin r + 0 O )dr. (3) 

J l Jo 

If the in-plane field is applied along the X-direction, the factor —h\\ sin r should be replaced 
by hypeos r (or, equivalently, one should operate the transformation p —> 1/p). From 
Eq.(j3]) with e 2 = 1 (p = 0) we easily recover the Fraunhofer-like magnetic diffraction pattern 
(MDP), | sin h\\/h\\ \ = |Sine /in |, typical of small linear junctions. 

The critical current, / c , of a quarter-elliptic junction can be found by maximizing 1 — Eq.([3]) 
with respect to 0 q: 


I c (±h \|) = 





9. 



r^t 2 


I"k/2 


/ X(r) sin(/iii sinrjdr 

+ 

/ X(r) cos(/i|i sinr)dr 


Jo 


Jo 


( 4 ) 


where / c (0) = J C WL is the junction zero-field critical current and E(e 2 ) = E(7r/2,e 2 ) is the 
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complete elliptic integral of the second kind of real argument e 2 < 1. Eq.Q is very general 
and, by properly adjusting the integration limits, applies to a junction shaped as any arc of 
ellipse with r in an arbitrary interval ['r 1 ,r 2 ]. 


B. Threshold curves computations 


Unfortunately, no analytical closed form exists for the definite integrals in Eq.Q which 
means that, in general, the MDP of a quarter-elliptic junction has to be computed numer¬ 
ically. The theoretical magnetic dependence of the normalized critical current, i c (±h n) = 
J c (±h|j)// c (0), of a small quarter-elliptic annular junction is shown for different p values 
in Figure [2} We observe that for quarter-circular junctions (p = 1), Eq.Q results in a 
Fresnel (or near-held) MDP typical of short Josephson junctions with asymmetric boundary 
conditiond 16 ^ 




In the above expression J n and H n are, respectively, the n-th order Bessel function of the 
Erst kind and the n-th order Struve function. We have also introduced the Fresnel’s Sine 
and Cosine integrals defined, respectively, by: S ( a / 2 / i ||/ 7 t ) = y/2h\\/n sin h\\x 2 dx and 
C(y/ 2h \\M = y/2h\\/n fo coshjjX 2 dx (h\\ > 0). In the limit h\\ —>■ 0, S ~ 0 and C ~ 
a/2/i||/7t, so that, as required, i c (0) = 1. In the opposite limit, i.e., for h\\ —* oo, S ~ C ~ 1/2 
and i c (h\\) asymptotically decreases as 1/ It is also interesting to observe that for a very 
prolate annulus, b»a , i.e., —e 2 ~ p 2 >> 1, the normalize MDP reaches a p-independent 
shape: 


i c (±h||) = + Hi 1 (/i||), 


which for large fields decreases as y/ir/h\\. 

For a semi-elliptic annular junction with r G [— 7t/2, 7t/2], taking into account the (j) — (j)o 
symmetry, Eq.Q reduces to: 


i c (±h\\) = 


E(e 2 


r ir/2 

X(r) cos(/i|| sin t)cLt 


(5) 


This expression also applies to full ellipses for which the parameter r spans over a 2tt range 
(as far as we disregard trapped Josephson vortices). The MDPs in Eq.([5]) are shown in 
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FIG. 2. Magnetic diffraction patterns of quarter-elliptic junctions for different values of the axes 
ratio, p. 

Figure [dj a) for different p values. For linear junctions (e 2 = 1), we again recover the 
expected Fraunhofer-like MDP, |Sinch||| = |jo(fy|)|, where j 0 is the zero order spherical 
Bessel function of first kind. For e 2 = 0 we end up with the Bessel-like dependence |J 0 (/in)| 
of circular and semi-circular junctions which first nulls for h\\ — 2.408 ~ n / {l-\-n~ l ). At last, 
for a very prolate ellipse, b»a , the normalize MDP reaches the limit p-independent shape 
7r|H_i(/i||)|/2. We observe that, as p increases, the secondary lobes grow, i.e., the MDPs 
look more and more like those of a linear junction with the Josephson current density, J c (r), 
peaked at the edges 15 (r = ±vr/2). 

It is evident that for an elliptic junction in a uniform held the minima in the magnetic 
pattern are not integer multiples of the first one, although they are (almost) equally spaced, 
the separation between two contiguous minima being about n. It is also worth to note that 
while the axes ratio p changes from 0 to oo, the first critical held, h?,, i.e., the hrst zero of the 
MDP only changes from n to 2 1 ' . This suggests that the dependence of /i| on p involves the 
complete elliptic integral of the second kind, E(fc 2 ), that spans the range [1,7r/2]; Figure^b) 
compares the numerically found values (dots) of hjj with the empirical expression (solid 
line) hjj(p) = 2E(ArcTanh 2 p) together with the percentage relative difference between the 
numerical findings and the proposed dependence (dashed line referred to the right scale). 
In passing, we remark that the integrals in Eq. (J5]) can be accurately reproduced by the 
following empirical functions: 

*c(=tfy|, p) — T(1 + a) ^ \ J A h \\)\’ ( 6 ) 
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(a) 


(b) 


FIG. 3. (a) Threshold curves of semi-elliptic and elliptic annular junctions according to Eq.|5]); 
n=0 indicates the absence of trapped Josephson vortices.(b) The normalized critical held /tjj vs. p: 
the dots are the numerical values, while the solid line is the function 2E(ArcTanh 2 p)\ the dashed 
line shows the relative error (right scale). 


where T is the Gamma function of real argument and 2cr(p) = (1 + n 1 )hUp) — n. 


C. Periodic boundary conditions 

For any (real) values of a and 6, it is possible to find a value, c, and a number, u, such 
that a = ccosliz/ and b = esinh //; c = ±a/ a 2 — b 2 are the abscissae of the ellipse’s foci. 
For a circle a — b, then c = 0 and v = oo; for a segment, 6 = 0, then c = a and v = 0. 
Furthermore, when |6| > |a|, then the foci lie on the E-axis and v = ArcTanh p is a complex 
number. (iy r) forms the so-called planar elliptic coordinate system. In general, the periodic 
conditions for the two-dimensional field (ft around an EAJTJ are written as: 


(ft(u, t + 27r) = (ft[y , r) + 27m, 

(7a) 

0 t(r T + 27t) = (ft T (u , r), 

(7b) 


n being an integer number corresponding to the algebraic sum of Josephson vortices (or 
fluxons) trapped in the junction at the time of the normal-to-superconducting transition. 
Eqs.(7a) and (7b) were derived in Ref. 18 for ring-shaped junctions and state that observable 


quantities such as the Josephson current (through sin0) and the magnetic field (through 
(ft T ) must be single valued upon a round trip. The net number of trapped fluxons comes 
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FIG. 4. Magnetic diffraction patterns of elliptic annular junctions as given by Eq.(|9]) with n = 1 
for different values of the axes ratio, p (compare to Figure [3](a) for n = 0). 


out to be the algebraic difference between the number of flux quanta associated with the 
fluxoids in each electrode 24 ! (f) in Eq.([2j is 27T-periodic, therefore, to implement the periodic 
conditions, one only has to add a term, nr, that accounts for n distributed 27r-kinks: 


4>(t) = h\\ sin r + nr + <j) 0 , 


( 8 ) 


Being (j) still an odd function (disregarding </> 0 ), we can use again Eq.Q to derive the most 
general expression for the threshold curve of a short EAJTJ: 


i c (±h\i,n) = 


2E(e 2 


X(t) cos(/i|| sin r + nr)dr 


( 9 ) 


The MDPs of Eq.Q with n — 1 are plotted in Figure [5] for different values of p. As we 
progressively decrease the ellipse minor semi-axis b (while keeping the major semi-axis a 
constant) an EAJTJ tends to a 2W x 2 a rectangular junction with a slit in the middle; with 


n 


= e 2 = 1, Eq.(9) reduces to: i c (±/i||,l) = 


sin/i||//i||— cos/iy 

/in 


= where j\ is the order 


one spherical Bessel function of first kind. In general, if n Josephson vortices are trapped in 
a slit realized in a linear junction, then i c (±h\\,n) = |j n (/i|i)|, in analogy with ring-shaped 
junctions for which i c (±h\\,n) = |J n (/i||)|. Ultimately, generalizing Eq.([6|, we have proved 
that the integrals in Eq. (J9]) can be approximated by normalized functions which involve the 
Bessel function of nonintegral order: 


i c (±/i||,p, n) = T(1 +n + a) 



To correctly interpret Figures [2j |3](b) and |3| one should keep the semi-axis a constant and 
change p through the semi-axis b. If a changes, then the held normalization changes. 
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III. LONG ELLIPTIC ANNULAR JUNCTIONS 


In this section we derive the appropriate partial differential equation (PDE) for an EAJTJ 
in an external magnetic field. The total tunnel current density is given by: 


J z = J c sin 6 + 


2 nR 


»t, 


where the last term takes into account the quasi-particle tunnel current assumed to be ohmic, 
i.e., R is the voltage independent quasi-particle resistance per unit area. The subscripts on 
(f) denote partial derivatives. Following Refs. 1 ^, a one-dimensional planar curved Josephson 
tunnel junction of constant width in the presence of a barrier-parallel external magnetic 
held, H e , is described by the following partial differential equation for </>: 


1 , • , / \ i 


ds 2 


uj 2 Rc s 


\ 2 j€ r € 0 R. 


+ 


A d H e ■ 0 


d 2 (pt 


ds 2 J r ds 


where oj 2 = 2ttJ c /& 0 c s , and c s being the specific junction capacitance. q(s) is the local 
normalized bias current density and A is the coupling between the external held and the 
held in the junction 118 . Here s is a curvilinear coordinate along the junction. It is well known 
that A j gives a measure of the distance over which significant spatial variations of the phase 
occur; the plasma frequency, uj p /2tt, is the oscillation frequency of small amplitude waves. 
v = ( x 0 /a 2 N,y 0 /b 2 N ) is the (outward) normal unit vector to the ellipse x 2 /a 2 + y 2 /b 2 = 1 at 
the point (a;o = a sin r, y 0 = b cost ) with N = \Jx\jaJ + y\ /ft 4 . Further, we can introduce 
the parameter c = uj p Xj which gives the velocity of electromagnetic waves in the junction 
and is called the Swihart velocity. The third term in the right-hand side takes into account 
the effect of the surface currents in the London layers, i.e., R s is the voltage independent 
surface resistance per area. Introducing the dimensionless loss coefficients cU 1 = c o p Rc s and 
fd = e r 6oU} p R s and normalizing the time t to uJp 1 , the last equation takes the form: 


,d 2 (p 


d 2 4>t 


~ ~ sin0 = 7(s) + afa - X + 


A dH„ 


( 10 ) 


3 ds 2 ru ' ^ Iy ' ' ' ri J ^ ds 2 J c ds 

where H u = H e ■ u is the component of the externally applied in-plane magnetic held 
normal to the junction perimeter. Therefore, in the experiments the magnetic held can be 


substituted by a properly chosen additional bias current 7 i(s) and vice versa. Eq.(10) is 
called Perturbed sine-Gordon Equation (PSGE). Because of its local form, it is quite general 
and holds for junctions of any geometrical shape. We note that the hrst and last terms on 


the right-hand side of Eq.(10), from a mathematical point of view, play the same role. 
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FIG. 5. Plot of p 2 Z 2 (p, r) sinr for 0 < r < tt and different p values. 

A. The sine-Gordon modeling 


We now want first to compute the normal component, H u , of a uniform magnetic field 
H = (0,i?||). H u = HyUy = /h|jX 1 (t) cos t is maximum at the ellipse poles, H u ( 0) = 
~H V (7r) = H\\, and vanishes at the ellipse equatorial points, r = ±7 t/ 2. Next, recalling that 
d _ drrd_ _ id ^ ie directional derivative of the normal field is: 

as as dr al dr 1 


dH u _ H\\p 2 


sm r. 


ds aZ A 

It can be shown that the second directional derivative is: 


( 11 ) 


ds 2 


a 2 X 2 


dr 2 


(1 -P 2 ) • d 

——— sm r cos r— 
X- dr 


( 12 ) 


Inserting Eqs.(ll) and (12) in Eq.(10), we end up with the following PSGE for an EAJTJ: 

„2 


+ 


P 


aI(p,T) 
A .7 


sin r cos t <p T + f3 ( 4> TTt + 


P 


X 2 


S111 T COS T (p Tt = 


1 2 


[4>tt + sin 4> + 7 (r) + afa] + h\\ A 


P* 


X\piT\ 


sm r, 


(13) 


where again h\\ = kclH\\. Eq.(13) states that for an elliptic annular junction the magnetic 
held enters directly into the PSGE in contrast to the case of linear junctions for which it 
appears only in the boundary conditions. Further, the different sections of the annulus feel 


different fields; diametrically opposed points feel opposite fields and the held term in Eq.(13) 
is out of phase with respect to the actual normal held. Moreover, the effect of a given held 
depends quadratically on the axes ratio p; the odd expression p 2 X~ 2 sinr, whose rms value 
is vX2. is plotted in Figure [H] for different p values. For p = X 2 (p, r) = 1, Eq.([l3| reduces 
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to well studied PSGE for ring-shaped junctions 118 . It is also worth to notice that, in the 


limit p —> 0, Eq.(13) reproduces the classical PDE for a linear junction. In fact, observing 


that Z 2 (0, r) = cos 2 r and that, as p gets smaller and smaller, p 2 Z 2 (r) sin r approaches the 
unitary impulse function, <5i(±7t/2) = ±1 and zero elsewhere (see Figure [5]), we have: 


4>rr + tan r 0 T = J cos r (0tt + sin 0 + 7 + a0t) - f3 (0 Trt + tan r 0 rt ) + hy Ahi(r). 

Now recalling that x = asinr, it is not difficult to derive that for any composite function 
0(x(r)) it is 0 TT + tanr 0 T = a 2 cos 2 r (j) xx \ introducing the normalized spatial coordinate 
£ = x/Xj, we recognize the well known PSGE for a linear overlap JTJ: 

- (frtt - sin 0 = 7 + acj) t - f3(f)^ t + h\\A5i(a/\j), 


with boundary conditions cj)^(a/Xj) = —a/Aj) oc h\\. In concluding, Eq.(13) is more 
general that it appears at a first glance. 


B. The static numerical simulations 


In what follows, we are interested in the static, i.e., time-independent solutions of Eq.(13) 


and, in order to have zero fluxons trapped in the junction, the conditions on the phase 
periodicity are: 0(r + 27r) = 0(r) and 0 r (r + 27r) = 0 r (r). The direct numerical integration 


of Eq. (13) with a = /3 = 0 poses large problems of stability due to the fact that there 


is no loss in the systemP^, to avoid this problems, we set a = 3 in order to have a fast 
decay towards a static solution (in real junctions a < 0.1). The term containing the surface 
loss was simply dropped to save computer time (/3 — 0). For the sake of simplicity, in our 
theoretical investigation the bias current was supposed to be uniform, y(r) = 7 , and A 
was set to 1 ; in real devices, depending on the specific electrodes geometry, non uniform 
bias distribution^H, self-held 22 23 and focusing 24 25 effects should be taken into account. In 
addition, we introduced the held normalization usual for long linear JTJs, /i( = H\\/ J c \j = 
h\\\jdj/ad m ] quite often in the literature the factor dj/d m > 1 has been ignoreef 26 . With 
symmetric bias and held terms, in the steady state regime with n = 0, the Josephson phase 


has the symmetry </>(± 7 t /2 + r) = </>(± 7 t /2 — r), therefore Eq.(13 can be integrated in the 


range [— 7 t/ 2 , 7 t/ 2 ] with boundary conditions 0 r (± 7 r/ 2 ) = 0. The commercial Unite element 
simulation package COMSOL MULTIPHYSICS (www.comsol.com) was used to numerically 
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FIG. 6. Numerical magnetic diffraction patterns of long EAJTJs for different values of the axis 
ratio, p, and two normalized perimeters: (a) t = 2n and (b) t = 47r. The magnetic held is 
normalized to J c Xj. 


solve Eq.(13) for different values of the normalized semi-axis, a/Xj , and of the mean axes 
ratio, p. Specifically, we have numerically computed the maximum value, i c , of the zero- 
voltage bias vs. /ij|. To begin with, we have checked that for a < O.lAj we were able to 
accurately reproduce all plots in Figure [3|a). 


For long EAJTJs we wish to confront the MDPs of samples having the same normalized 
perimeter, i = P/Xj = AaE(e 2 )/Xj, and different axes ratios; accordingly, for given d and 
p, we must set a/Xj = 1/AE{1 — p 2 ). In order to compare the numerical findings with 
the experimental data, we will limit our interest to p = 0.5 [E(0.75) ~ 1.21] and p = 2 
[E(—3) ~ 2.42], while the cases of well understood linear geometry, i.e., p = 0 ]E(1) = 1] and 
circular geometry, i.e., p = 1 [E(0) = 7r/2], are taken for reference. Pronounced deviations 
from the theoretical behavior of short junctions were found for l = 27T, as can be seen 
in Figure |6j( a) where the i c (h y) dependence is reported. For small fields a long linear JTJ 
behaves as a perfect diamagnet by establishing circulating screening currents which maintain 
the interior held at zero. This ’’Meissner” effect is not anymore possible in long curved 
junctions where the magnetic held does penetrate the barrier even when small; nevertheless, 
still the critical current decreases linearly with the applied held. We also observe that some 
ranges of magnetic held develop in correspondence of the pattern minima in which i c may 
assume two different values. In a fashion which closely resembles the behavior of long linear 
junctions, these values correspond to different configurations of the phase inside the barrier. 
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In fact, each lobe in Figures [6] is associated with a given vortex structure; at the very end 
of the first lobe a fluxon-antifluxon pair is present along the junction, the fluxon facing 
the antifluxon on diametrically opposed potential wells created by the magnetic held at 
the ellipse’s poles. In the successive lobes the magnetic held penetrates in the barrier and 
more pairs are nucleated at the ellipse’s equatorial points, in a way which closely recalls the 
behavior of the type II superconductors, even though the vortices we are dealing with are 
quite different from the Abrikosov vortices having a normal core. In the second lobe, for 
example, we start from a phase conhguration very similar to that at the right side of the 
hrst lobe in which one huxon-antihuxon pair is present in the barrier and we end up with 
two huxon-antihuxon pairs, the two bunched huxons facing the two bunched antihuxons 
on diametrically opposed potential wells. In order to trace the different lobes of i c vs hj,, 
it is crucial to start the numerical integration with a proper initial phase prohle. As / is 
increased, each higher lobe broadens further and more lobes overlap. This behavior is shown 
for £ = 47r in Figure (6|b) where also the third lobe was computed. 

Let us recall that in Figures j6ji c is plotted versus h'^ (rather than hy) which corresponds 
to a normalization of the magnetic held to J c Aj, as it is usual for long linear Josephson 
junctions for which, as £ increases the critical held saturatej F( 18 l 24 l 2 § E9] a t 2J c Xj. For one¬ 
dimensional ring-shaped junctions it was found 13 3 ® that the critical held is proportional to 
the ring radius, r, as far as r » Xj. Our numerical simulations allowed us to reach the 
conclusion that for long elliptic annuli the critical held increases as ( a/Xj) p . 


C. The self-field effects 

The analysis of long EAJTJs would not be satisfactory if we neglected the effects of the 
magnetic field generated by the junction feed current; for the sake of generality, we will 
treat them only on a qualitative basis, since this effects drastically depend on geometrical 
details. It will turn out that the following considerations are very useful to interpret the 
experimental finding reported in Section V. Let us consider an EAJTJ fed by a d.c. bias 
current, /, smaller that the critical current, / c , so that it is in the stationary (zero voltage) 
state. As depicted in Figure [7|a), / is applied in the positive X-direction: it enters the 
junction, say, from the wiring electrode on the left side, gradually splits in the two arms of 
each electrode, recombines on the junction right side and leaves through the bottom electrode 
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FIG. 7. (a) Schematic representation of the in-plane normal magnetic self-held generated by the 
tangential supercurrent, I t , flowing in the top electrode of a biased Lyngby-type EAJTJ with 
H = 0. (b) Vertical cross section (not to scale) of our Lyngby-type elliptic annular Josephson 
junction. 


(on the right). It is convenient to focus on one electrode at a time and we consider first the 
top/wiring electrode. In absence of an external magnetic held, the symmetry of the system 
requires that the supercurrent, 7f(r), flowing along the top electrode within a distance of A t 
from the Josephson barrier is null at the equatorial points and equally distributes in each 
annulus branch; this supercurrent increases as we move towards the poles where it must be 
equal to J/4, since I = A(0) + A(0) + A( n) + A( tt) and A(0) = I t (0) = A( 7 r) = A( 7r). It 
is important to stress that It hows within a width, W t , that is constant all along the ellipse 
perimeter. This can be better understood looking at Figure [7](b) which sketches the vertical 
cross section along the equatorial line of our Lyngby-type EAJTJs. The oblique arrows in 
Figure [T)(a) indicate the strength and direction of the normal component of magnetic held 
associated (and proportional) to the tangential supercurrent I t that is inversely proportional 
to Wt- We can make a similar reasoning for the current, A> flowing in the bottom electrode 
(within a depth A&) with the caveats that A hows on the opposite side of the barrier so 
its normal held has opposite direction and, above all, A hows within a width equal to (and 
where possible larger than) IJJ that, in turn, is larger than W t . It means that, on the 
average, the normal held induced by A is much smaller that the one generated by A (this 
would still be true even if W t = W 6 ). In our samples W t = 0.3UJ,; therefore, for a qualitative 
understanding, we will only consider the self-held contribution due to A- From Figure [7](a) 
it is evident that such normal held has a net V-component, while, in the average, the X- 
component is null. Further, the Y -component of the normal self-held is larger for oblate 
ellipses (p < 1) and vice versa. The V-component affects the prohle of the Josephson phase, 
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therefore, as / is increased we will reach a value that makes the junction to prematurely 
switch to the voltage state (the largest critical current always corresponds to the uniform 
phase distribution, 0(r) = sin _1 /// c ). In the presence of an externally applied magnetic 
held, the system symmetry is broken and, in general, the bias current splits in unequal 
parts giving rise also to a net X-component of the normal self-held which is larger for 
oblate geometries (p > 1). In a hrst approximation, as already heuristically suggested for 
ring-shaped junction^, the self-held effects in long EAJTJs can be taken into account in 


Eq.(13) by adding extra terms which simulate two uniform helds, one parallel, oc pycosr, 
and the other orthogonal, oc p^ysinr, to the bias current how. 

We should also add that, even neglecting the self-held effects, a bias current density which 
is uniform in r can never be achieved in AJTJs. If the bias current, I, were uniformly 
distributed along the height of the current carrying electrode, 7 (s(r)) would be peaked at 
the equatorial points of the ellipse. However, the current in a superconducting hat him 
mainly hows at its boundaries, therefore, more realistically, the 7 prohle is depressed at 
the equatorial points and is largest at the poles. Although symmetric, a non-uniform current 
density prohle, together with the self-fields, reduces the largest possible value of the critical 
current. I 11 the experiments, the ratio of the zero-held critical current, I C) 0 , to the current 
jump, A I g , at the gap voltage is a direct measure of this nonuniformity: the lower this ratio, 
the larger is the nonuniformity. 


IV. TRANSVERSE MAGNETIC FIELD H± 

An alternative way to modulate the critical current of a planar Josephson tunnel junction 
is to apply a magnetic held. Hi = (0, 0, H±), perpendicular to the junction plane®®!, which 
induces shielding currents in its electrodes. In turn, the shielding currents generate a local 
magnetic held whose normal component thread the Josephson barrier. The modulation 
amplitude drastically depends on the geometry of the electrodes and on how close to the 
barrier the shielding currents circulate. For rectangular junctions these effects have been 
already investigated both theoretically and experimentally 14 -^, for ring shaped junctions 
only magnetostatic simulations and experimental data exist, and a theoretical understanding 
is still lacking. Here, in addition, we will analyse how a transverse held acts on an EAJTJ; 
more specifically, we want to derive the normal component of the magnetic held induced by 
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(a) (b) (c) 

FIG. 8. Sketches of the bottom electrode of elliptic annular junctions having different geometries: 
(a) island geometry, (b) modified Lyngby-type geometry and (c) Lyngby-type geometry. The 
top/wiring electrodes are specularly symmetric. 

the the circulating shielding currents. We will consider three different geometries formed by 
specularly symmetric electrodes whose bottom electrodes are shown in Figures [8]( a)-(c). We 
will take in consideration that for window-type junctions the elliptic loops in the bottom 
and top/wiring electrodes have the same mean axes, i.e., the same perimeter and cllipticity, 
but different widths (respectively, W b and W t ). For the sake of simplicity, the analysis will 
be carried out assuming that no flux is trapped in the electrode loops. 



A. Island geometry 

The island or ^-biased 20 geometry, in which each electrode is fed by a very tin thread, 
is the easiest to analyze since in the Josephson region the electrodes have constant widths. 
Let us consider first the bottom electrode sketched in Figure |8](a). When a magnetic field, 
H l, is applied perpendicular to the loop plane, then a shielding current, I ar .t, = /iq H±A b /Lb, 
circulates in the loop to restore the initial flux 3 ®, where A b is the effective flux capture area of 
the loop in the bottom electrode and L b is the inductance of a thin-hlm elliptic (narrow) loop. 
What really matters is the magnetic flux <f>£, = /a 0 H±Ab applied to (although not threading) 
the loop; to a very high degree of approximation the capture area in this case is the loop 
inner area A b ~ nab. Applying the Ampere’s circuital law to the (tangential) shielding 
current, I cir ^, we easily derive the normal magnetic field, H ub = I cirb /2Wb = K b H±, where 
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Kb = HoAb/2WbLb is the conversion factor from the transverse to the barrier-parallel held 3 ' . 
If we now consider the top/wiring electrode, for symmetry reasons and considering that 
the shielding current circulates on the opposite side of the barrier, we have that H U)t = 
—K t H±. In the case of equal width annuli, then Kb = K t and the total normal held, 
H Uy b + H U:t , is null meaning that the held generated by the currents circulating in the top 
electrode fully compensates that produced by the shielding currents induced in the bottom 
one. However, in general, A K = Kb — K t ^ 0 and the total normal held is uniform 
and proportional to the transverse held, H u = AK H±, therefore, for short EAJTJs it 
gives rise to Fraunhofer-like threshold curves. To demonstrate this we resort to Eq.Q 
in elliptic coordinates, (<90/<9iy <90/<9r) = KaI(r)(iJ r , — H u ), which, in our one-dimensional 
approximation, yields d^/dr = —kclH u Z{t) = —kclAKHKZ{j). By integrating dcft/dr and 
introducing the dimensionless held hj_ = na,H±, we have: 

0(r) = -h ± AK E(t, e 2 ) + 0 O . (14) 


Reiterating the calculation of Sec.IIC, with 0 as in Eq.(14), we get the MDP of a short 
island-type EAJTJ in a transverse held: 


A(±h_ l) = 


E(e 2 


r 7r/2 


X(r) cos [h±AK E(r, e 2 )] dr 


E(e 2 ) 


fE(e 2 ) 

cos h i zdz 


'0 


sin h±AK E(e 2 


ftiAA'Efe 2 


= I Sine h±AK E(e 2 ) I . 


(15) 


Being its derivative is null, according to Eq.(10), a uniform normal held produces no effects 
on electrically long EAJTJs. 


B. Modified Lyngby-type geometry 

We now analyse the modified Lyngby-type geometry in which the left (right) side of the 
bottom (top) electrode is a semi-infinite plane; the bottom electrode of the configuration is 
shown in Figure |8](b). On the right side of the bottom electrode, the top/wiring electrode 
acts as a ground plane and squeezes the held lines generated by the screening currents in such 
a way that most of the magnetic energy is confined between the electrodes. The inductance 
per unit length of such superconducting microstrip transmission line 35 ® is Cb = poAfe/HJ,, 
where A b is the bottom current penetration! 23 39 ^ coth ^ + y tanh We claim that the 
contribution to the loop inductance of the left side of the bottom electrode is negligibly small, 
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therefore the loop inductance is given by the product of the ellipse mean semi-perimeter by 
the bottom inductance per unit length, L b = PC b /2. This inductance is much smaller than 
that of a free standing elliptic loop used for the island-geometrjJ 24 l 

Being W b » d m) for 0 < r < 7r, the normal magnetic held penetrating the barrier is 3 ' 10 
H V)b = I C ir,b/W b = K b H±, where now the held conversion factor 35 is K b = 2A b /PA b • We note 
that K b is independent of W b , since the circulating current linearly increases with the strip 
width. The shielding currents mostly how along the outer electrode boundary, therefore, 
for —Ti < r < 0, no normal held threads the Josephson barrier. If we now focus on the 
specularly symmetric top/wiring electrode, we have that H U)t = —K t H± for — tt < t < 0 
and zero elsewhere. Here K t = 2A t /PA t , with A t = A coth W _|_ A t a nh W = dj — A t . 
Summing up the contributions from both electrodes, the total normal held is given by a 
27r-periodic step function having different amplitudes in the two half-periods: 


H , 


—K f if — 7T < r < 0 
K b if 0 < r < 7T. 


(16) 


Starting from Eq.(16), it is not difficult to derive that for short elliptic junctions with the 


modified Lyngby geometry in a transverse magnetic held: 


i c (±h ± ) = - 


r,. ^d m ^ b 7Tcf m $j 

bme ——z-b Smc 


Afod'o At<h 0 

which, in the case of equal huxes, again reduces to a pure Fraunhofer pattern that only 
depends on the coupled hux (and not at all on the junction cllipticity). It is worth to 
point out that, due to hux focusing effects, the effective capture area of this geometry is 
much larger than the ellipse inner area. As far as long junctions concerns, the normal 


held discontinuities at the ellipse poles result in two Dirac terms in the PSGE in Eq.(10), 
corresponding 20 to local extra conditions on the spatial phase derivative. 


C. Lyngby-type geometry 

Let us consider the bottom electrode of a Lyngby-type EAJTJ sketched in Figure [8](c). 
The loop width is constant and equal to W b , on the right side and increases when we move 
to the left side. With respect to the previous geometry, this makes the transition to zero 
of the bottom normal held smoother as we cross the ellipse poles. The calculation of how 
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FIG. 9. Plot of /(p, t) = 1 + pZ 4 (p, t) sin r for — 7 r < r < 0 and different p values. 


this field vanishes when the circulating current moves away from the barrier is not an easy 
task, but, clearly this effect occurs faster for prolate junctions. Being W b « b, we will 
suppose that the held decay only depends on p and can be qualitatively described by the 
function /(p, r) = 1 + pl _1 (r) sin r plotted in Figure [ 9 ] for different values of p (in the range 
— 7 r < r < 0). Summarizing, the normal held, //^(r), generated by the screening currents 
circulating in the bottom electrode is given by the continuous expression: 


H Vtb (r) 

H 1 


K b f(p, r) if - 7r < r < 0 
K b if 0 < r < 7T. 


(U) 


Analogously, for the specularly symmetric top electrode it is: 


H v , t {r) 

H 1 


—K t if — 7T < r < 0 

-K t f(p, -r) if 0 < r < 7T. 


(18) 


Whatever are the electrode inductance, the conversion factors are 42 K b) t = poAb } t/2Wb,tLb,t- 
The total normal held is H v = H U)b + H U)t , and is given by the continuous function: 


H A r ) _ J T )~ K t if - 7T < r < 0 

K h - K t f(p, -t ) if 0 < r < 7T. 


H 1 


(19) 


In the case of symmetric electrodes (K b — K t = K), the above expression greatly simplifies 
to H v (t) = H^KpZ~ l {j) sinr which is the normal component of a uniform in-plane held 
KH± applied in the direction of the Waxes, which is the direction of the current carrying 
electrodes. Remarkably, the critical current of both short and long Lyngby-type (symmetric) 
EAJTJs in a transverse held is expected to modulate exactly as described in Sections II and 
III for an in-plane magnetic held. This result is supported by magnetostatic simulation^ 14 
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FIG. 10. Sketch of an elliptic annular Josephson junction with ratchet geometry. The bottom 
electrode is in gray while the top/wiring electrode is black. 


showing that for Lyngby-type ring-shaped junctions in presence of a transverse field the 
normal magnetic held has a sinusoidal dependence on the polar angle 9, independent of 
the annulus radius. Later on, it was experimentally proved 13 (although at that time not 
explained) that the held conversion factor K of Lyngby-type ring-shaped junctions increases 
linearly with the mean radius, r; this is consistent with the capture area proportional to 
r 2 and an inductance proportional to r (as in the case of the modified Lyngby geometry). 


Further support of the validity of Eq.(19) will be provided by the experimental data reported 
in Section V. 


With Kf, ^ K t , Eq.(19) results in an asymmetric held prohle, H v {—t) ^ —H u (t). The 
consequences resulting from the asymmetric boundary conditions imposed by a non-uniform 
external magnetic held at the extremities of both short and long linear Josephson junctions 
have recently been investigated; the held asymmetry is responsible for a degeneracy of the 
(extrapolated) critical held, Hp that was numerically demonstrated^ 6 and experimentally 
verified® Asymmetric held prohles in long Josephson junctions were recently looked foi 19 ® 
in order to exploit the rectifying property® of a ratchet potential. It is easy to recognize that 


the configuration depicted in Figure 10 (bottom electrode in gray and top/wiring electrode 
in black) implements the ideal step-like deterministic ratchet potential with the normal held 
being constant on the left side (K t H±) and null on the right side. 
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FIG. 11. (Color online)(a) Optical image of one of our samples; (b) Exploded sketch of our layout; 
for clarity, the top/wiring layer (in blue) is shifted slightly upward with respect to the bottom layer 
(in red). 


V. THE MEASUREMENTS 


A. The samples and the experimental setup 

Using the well known and reliable selective niobium anodization procesd^ 5 we have realized 
high-quality window-type Nb/Al-AlOx/Nb EAJTJs. The details of the trilayer deposition 
and of the fabrication process can be found elsewhere®. As shown by the photograph 
in Figure [TT|(a.) , two Lyngby-type elliptic annular junctions, named A and B, were made 
having the same eccentricity, but rotated 90° relative to one another. The junction minor 
and major inner axes were, respectively, 140 and 280 /im and the nominal annulus width 
was W = 10 /im\ however, due to different anodization conditions from wafer to wafer, 
the effective width was in the range of 7-9 nm. Figure [TT|(b) shows in more details the 
chip layout in which the elliptic tunnel barrier is sandwiched between two simply connected 
superconducting layers; the top/wiring layer (in blue) is shifted slightly upward with respect 
to the bottom layer (in red). 

Niobium anodic oxide (80 nm thick) and an extra dielectric layer made of rf-sputtered silicon 
dioxide (120 nm thick) provided the electrical insulation between the base electrode and the 
wiring film around the junction area. This 200 nm thick insulating by-layer outside the 
elliptic ring, also called idle region, alters both the static and dynamic properties of the 
junctions; on linear one-dimensional long junctions it has been proved by both numerical 
simulation 47 and experiments® that, as far as the static properties are concerned, the only 
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effect of the idle region is to increase the magnetic energy stored in the fluxons, i.e., it 
introduces a scaling factor on the field strength. The different widths of the narrow part of 
the base and top electrodes, respectively, W], = 20 and W t = W = 6 /im, result in a lack of 
full specular symmetry of our samples. The different thicknesses of the base and top/wiring 
layer also contribute to the system asymmetry. For our samples the thicknesses of the base 
and top/wiring Nb electrodes were, respectively, db = 200 nm and d t = (40 + 500) nm ~ 
6 A 7 V 6 , with A m{T = 4.2/1) = 90 nm, resulting in magnetic and current penetration depths 
of, respectively, d m ~ 162nm and dj ~ 182 nm ~ 1.12 d m . Furthermore, the bottom and top 
current penetration depths were, respectively, A?, « lOOnrn and A t = dj — A^ ~ 82 nm ~ 
0.82A;,. The critical current density of our samples was measured on electrically small cross- 
type junctions realized during the same deposition batch on different chips; at T = 4.2 K, 
we found J c = 53 A/cm 2 corresponding to A j ~ 52 pm. Figure [T)(b) also shows that in our 
samples there is a 5 /im wide idle region only on the outer left side of the junction; taking 
into account this asymmetric idle region, it is A j ~ 58 /im on the junction left side. For our 
calculation we will use the average value A j ~ 55pm. 

Our setup consisted of a cryoprobe inserted vertically in a commercial LHe dewar. The 
cryoprobe was magnetically shielded by means of two concentric Pb cans and a cryoperm 
one; in addition, the measurements were carried out in an rf-shiclded room. The external 
magnetic field could be applied both in the chip plane or in the orthogonal direction. The 
chip was positioned in the center of a long superconducting cylindrical solenoid whose axis 
was along the Y -direction (see Figure [I]) to provide an in-plane magnetic field, i/||. The 
transverse magnetic field, H±, was applied by means of a superconducting cylindrical coil 
with its axis oriented along the Z-direction. All measurements were carried out at T = 4.2K. 


B. In-plane magnetic diffraction patterns 

In this section we report the measurements carried out on junctions A and B of Figure [IT] 
having, respectively, Pa = 2 and ps = 0.5, according to our notations (see Figure [I]). For 
both junctions the mean perimeter is P = 4E(1 — p 2 B ) as ~ 4 x 1.211 x 145 pm ~ 700 pm, 
i.e., much longer than the Josephson penetration depth, £ — P/Xj ~ 13 (it is E(l — p 2 B ) = 
PbE(1—Pa) ~ 1-21). A large number of such samples were investigated and they all showed 
not only the zero-field critical current, J Ci0 , but also the maximum critical current, I c ,maxi 
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FIG. 12. Experimental magnetic diffraction patterns of our EAJTJs with an in-plane magnetic 
field: (a) junction A and (b) junction B. The extrapolated dotted lines help to locate the critical 
fields. 


considerably smaller than about the 70% of the current jump at the gap voltage, AI g , 
typical of short Nb / Al-AlOx-Al / Nb junctions. As anticipated in Section IIIC, this is the 
first signature of a non-uniform bias current distribution and, more importantly, of the 
self-field effects. 

Figures [l2|(a)-|T2)(b) display the MDPs, respectively, of samples A (A I g = 5.7 mA) and B 
(AI g = 5.5 mA) with an in-plane magnetic field applied in the direction perpendicular to the 
bias current flow. At a first glance, we observe that both datasets are quite different from 
what is expected for a long EAJTJ with a similar normalized length (see Figure [6](b) ). With 
real devices, the measurements of maximum supercurrent as a function of the external field 
often yield the envelop of the lobes, i.e., the current distribution switches automatically to the 
mode which for a given field carries the largest supercurrent. Sometimes, for a given applied 
field, multiple solutions are observed by sweeping the junction current-voltage characteristic 
many times. Nevertheless, the (first) critical field T/| can still be obtained extrapolating 
to zero the MDP first lobe (see dotted line in Figures [T2)(a)-|T2|(b)). Furthermore, while all 
the theoretical threshold curves derived in the previous sections were implicitly meant to be 
symmetric with respect to the bias current and/or magnetic field inversion, the experimental 
MDPs only retain the symmetry with respect to the simultaneous inversion of the bias 
current and of the magnetic field. The pattern skewness is mainly to be ascribed to the Y- 
component of the normal self-field; in the second and fourth MDP quadrants this component 
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(a) (b) 

FIG. 13. Experimental threshold curves of our EAJTJs in a transverse magnetic field: (a) junction 
A and (b) junction B. The dotted lines help to locate the critical fields. 


adds to the external field, while in the first and third quadrants it partially compensates 
the applied field; this results in different slopes of the first lobe and therefore in different 
(extrapolated) critical fields. Our data are consistent with the postulate in Section IIIC, of 
the current-perpendicular self-held being larger for prolate EAJTJs. In different words, the 
in-plane MDP of our samples would not have been skewed, if the magnetic held were been 
applied in the direction of the bias current. Furthermore, from the numerical simulations 
of Figure |6](b) we expected the critical helds of junctions A and B to be, respectively, 
~ 5J c Aj ~ 145 A/m and ~ 2.4J c Aj ~ 70 A/m which are close to the measured average values 
(respectively, 130 and 65 A/m). The larger I c , max /AI g ratio observed for junction B, can be 
explained by smaller current-parallel self-fields in oblate rather than prolate EAJTJs, but 
is also consistent with the supposition that, as the height of the current carrying electrodes 
gets smaller, the current density prohle becomes more uniform. 


C. Transverse magnetic diffraction patterns 

Figures [I3fa)-[l3l(b) display the I c vs. H i dependences of, respectively, the same junctions 
A and B reported in Figures |l2)(a)-Jl2|(b). For sample A we observe that the transverse MDP 
is practically indistinguishable from its in-plane counterpart (apart from a held factor scale); 
this confirms our finding of Section IV that for Lyngby-type annuli a transverse held is 
equivalent to an-in plane held applied in the direction of the current how and at same time 
substantiate the presence of a current-parallel self-held (otherwise the transverse pattern 


25 



































would not be skewed). The comparison of the transverse and in-plane MDPs of sample 
B further supports this picture, if we recall that the effect of the current-parallel self-field 
is reduced, if not negligible, in oblate junctions. Now the transverse critical field is larger 
for junction B (7.5 A/m), rather than junction A (4.1 A/m)] this was expected considering 
that a 90° rotation of the in-plane field corresponds to a transformation p —)■ 1/p. In other 
words, in presence of a transverse magnetic field, the system symmetry is broken along the 
direction of the current flow. This makes sample A about 32 times more sensitive to Hj_ 
than to H \\, the field conversion factor being about 8.6 for junction B. 


D. Vortex trapping 

It was possible to trap Josephson static vortices (fluxons) on a statistical basis by means 
of fast coolings of the samples through their superconducting transition. The trapping 
probability is know to grow with the speed of the normal-to-superconducting transition 49 50 . 
After a successful trapping procedure the zero-voltage critical current is considerably smaller 
and a stable finite-voltage current branch, called zero-field step, appears in the junction 
current-voltage characteristic indicating that the bias current forces a single fluxon to travel 
along the ellipse perimeter in the absence of collisions. The corresponding flux quantum must 
be trapped in the superconducting loop formed by either the bottom or top electrode. For our 
samples the current branch associated with one fluxon had an amplitude larger than 1 mA 
and an asymptotic voltage V± ~ 29 pV which results in an average speed, PV l/$o ~ 10 7 m/s, 
considerably smaller than the Swihart velocity, 1.5 x 10' m/s, typical of all-Vfe junctions 
evidencing, once again, indicating that the fluxon travels in the periodic potential^ generated 
by the bias current. The dynamic properties of EAJTJs involves the nucleation, propagation 
and interaction of more than one fluxon or fluxon-antifluxon pair and will be discussed in a 
future work. 


VI. CONCLUSIONS 

The static properties of one-dimensional elliptic annular Josephson tunnel junctions have 
been investigated theoretically and experimentally. Both short and long junctions were 
considered in the presence of an in-plane as well as a transverse magnetic field. For short 
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annular junctions we derived and computed the dependence of the critical current on a 
uniform in-plane magnetic held; it is found that the I c vs. H dependence is determined 
by the sample cllipticity, although the first critical held can only change by 50%. Different 
geometrical conhgurations have been considered for studying the effect of a magnetic held 
applied perpendicular to the ellipse plane; we established that for a Lyngby-type EAJTJ a 
transverse held emulates an in-plane held applied along the direction of the current how. 
Further, we derived the proper perturbed sine-Gordon equation to describe both the statics 
and the dynamics of the phase difference across the barrier of an EAJTJ. The static solutions 
of the partial differential equation were numerically computed and we found that, for a 
given held, different phase prohles are possible depending on the number of static huxon- 
antihuxon pairs nucleated at the ellipse poles where the derivative of the normal held is 
largest. We also evaluated the static self-held effect in long EAJTJs. Two planar structures 
characterized by not simply-connected electrodes have been considered in the experiments. 
Experimental data on high-quality long Nb/Al-AlOx-Al/Nb EAJTs basically confirm the 
numerical predictions, provided the effects of the current limiting static self-held are taken 
into account. A transverse magnetic held is demonstrated to be several times more efficient 
than an in-plane one to modulate the junction critical current. For a given inner area, 
prolate EAJTJs (for which the inversion of the periodic potential occurs faster) are more 
efficient for the in-plane to transverse held conversion. 


Unbiased elliptic annular junctions inherently have specular symmetry with respect to their 
principal axes: quite obviously an in-plane magnetic held breaks the system symmetry 
along its direction. In this paper we have demonstrated that a transverse held breaks the 
symmetry along the direction of the current carrying leads; furthermore, in long EAJTJs 
the bias current itself also generates non-symmetric conditions, the asymmetry being more 
pronounced along the current direction for oblate ellipses and vice versa. Among other 
things, we have also suggested a simple geometrical configuration in which the magnetic 
held coupled to the elliptic barrier lacks reflection symmetry, so accomplishing a nearly 
ideal rectifying potential in which a soliton is accelerated only in one half of the junction 
perimeter. The soliton dynamics in EAJTJs will be treated in a forthcoming paper. 
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